Topological properties of the Cosmic Microwave Background polarization map 
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Rapid progress has been made in observations of the temperature anisotropics of the Cosmic Mi- 
crowave Background (CMB). These observations have enabled cosmologists to characterize the state 
of the universe at recombination, and observational efforts are now being directed towards obtaining 
the polarization map of the CMB. Here we draw an analogy between the CMB polarization map 
and nematic liquid crystals, pointing out that the two have similar defects. Making use of known 
results in the theory of defect formation we predict the statistical distribution of defects in the CMB 
polarization map and provide a novel tool for probing the CMB at large angular scales. 



A polarization map of the CMB wiU consist of head- 
less vectors laid on the surface of last scattering. One 
can view the polarization map as rod-shaped molecules 
laid out on a two-dimensional sphere (5^), the sky. This 
picture immediately connects the CMB polarization map 
with nematics where the molecules are also rod-shaped. 
The topological properties of nematics have been stud- 
ied for many decades ^ . The very topological properties 
that are relevant to nematics will also apply to the CMB 
polarization map, as first examined in Ref. 0]- 

At first sight, it might appear that the analogy is 
doomed because neighboring molecules in nematics in- 
teract, while CMB polarization vectors do not. However, 
the processes at last scattering that produce the polar- 
ization are described by interacting field theory and so 
neighboring polarization vectors, just as neighboring liq- 
uid crystal molecules, tend to align. The CMB polariza- 
tion map is just like a snapshot of a thin film of nematic 
liquid crystal. 

The mathematical description of the CMB polarization 
map0 in a local patch of the sky is in terms of a 2 x 2 
symmetric traceless matrix which we denote by P: 



Q U 
U -Q 



(1) 



We define a normalized "order parameter" as: M = 
^P/ y^Q^+lP. The matrix M can also be written in 
terms of a unit vector n as: 



FIG. 1: Fundamental defects of charge ±f/2. Each dash 
respresents the linear polarization of the CMBR at that point. 
The X marks the position of the singularity where the linear 
polarization vanishes for topological reasons. 
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So the space of all matrices M is given by the space 
(S^) of all two-dimensional unit vectors but with the Z2 
identification under h — > — n. The normalized field M 
therefore defines a unit vector field with no arrows, also 
known as a "line field." 

Since M lies in /Z2, a space that has incontractable 
loops, the line field will have singularities, also called "de- 
fects." Some of these singularities are shown in Fig. 1. 
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FIG. 2: Defects of charge ±1. These can be constructed by 
combining two fundamental defects shown in Fig. 1. 



The radial and tangential singularities commonly dis- 
cussed in the CMB literature and shown in Fig. 2 can 
be constructed by combining two of the fundamental sin- 
gularities. 

In Euclidean space, each singularity has a degree W 
which is given by: 



Wr^ — (f dx' Tr(Mea,M) 
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where T is a contour around the location of the singular- 
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ity, and 
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In terms of the functions Q and U , we simply have 

Wr^^i dx'd.a (5) 

where a = (l/2)tan~-'^(t//(5)- The fundamental defects 
shown in Fig. 1 have degrees ±1/2, while the defects 
in Fig. 2 have degrees ±1. In general we can take F 
to be any closed contour on which the polarization is 
everywhere non-vanishing, and then Eq. ^ tells us the 
number of singularities within that contour. (We assume 
that all the singularities are isolated.) 

Equation Q can also be used on the surface of a sphere 
provided we first parallel transport M along a geodesic 
from points on F to a point within F. In other words, to 
compute the degree, we have to parallel transport M to 
a common tangent space. 

The Poincare-Hopf theorem tells us that the alge- 
braic sum of degrees of all singularities, also called the 
"index" of a vector field, on a compact manifold is equal 
to the Euler number of the manifold. By explicit eval- 
uation we can show that the Poincarc-Hopf result also 
applies to a line field on . Hence the CMB polariza- 
tion map must necessarily have index 2 since the Euler 
number of 5^ is 2. 

The distribution of defects produced at a phase tran- 
sition has received much attention over the past few 
decades . A key result is that if we consider the length, 
L, of F to be much larger than the correlation scale, ^, 
of the matrix M, then we can write @ 

WT,rms = {W^Y'^ ^c[^^ (6) 

where c is a system-dependent constant but the exponent 
(3 is universal with the value (3 = 1/2. The reasoning 
leading to this scaling is very simple. Wy is the sum of 
random fluctuations in a and the number of elements in 
the sum goes like L/^. Therefore the root mean square 
value of Wt is proportional to \/~lJ\- 

Satellite observations of the CMB 0, show that the 
temperature anisotropy correlation function, C{9), falls 
off very steeply with increasing angular separation, and 
becomes very small above an angular scale 6 = x ^ ^ 
few degrees. If this lack of large-scale temperature corre- 
lations is indeed real, and not some misunderstood sys- 
tematic effect, it is natural to expect that one should find 
a corroborating lack of correlation for polarization on all 
angular scales larger than x- After all, thermal and po- 
larization fluctuations are expected to arise from similar 
processes during recombination. Hence, on angular scales 
larger than the matrix M would be essentially uncor- 
related. This tells us that the areal density of defects on 
the sky will be approximately one per x'^ ~ 10 square 



degrees. Further, we expect (3 = 1/2 on angular scales 
larger than x- K the data shows P ^ 1/2 or exponen- 
tial decay with L, then that would indicate the presence 
of large-scale correlations in the CMB polarization. This 
observation would be equivalent to defect formation with 
"bias" in condensed matter systems. 

If the statistical distribution of the polarization is 
Gaussian, then all correlation functions of the polariza- 
tion matrix P can be calculated in terms of the two 
point correlation function. However, correlators such as 
(VFp) that are central to our discussion, cannot be writ- 
ten in terms of two point correlators of P, since they are 
given by correlators of the normalized matrix M which 
involves dividing P by y/Cp + lP. Alternately, {W^} de- 
pends on the two point correlation function of the angle 
a = (l/2)tan~^(t//(5) and this is not given by the two 
point correlation functions of U and Q. Hence the infor- 
mation provided by the distribution of defects is comple- 
mentary to that provided by the two point correlators of 
the polarization matrix usually discussed. 

There is a global constraint on the CMB polarization 
map, namely that "B-modes" are absent or highly sup- 
pressed HjUilllilllH and the reader may worry that this 
would affect the distribution of defects. However, the B- 
mode relates the variation of the polarization intensity, 
\/Q^ + W^, to the distribution of the unit vector n. The 
topology, however, does not depend on the polarization 
intensity, only on the direction n. Hence for a fixed dis- 
tribution of n, we would need to find a/Q^ -I- C/^ such 
that the B-modes vanish. One can show that even pur e 
E- modes exhibit the defects shown in Figs. 1 and 2 jla |. 

As polarized CMB radiation propagates from the last 
scattering surface to us, intervening effects such as weak 
lensing and Faraday rotation can affect the polarization 
map. As long as these effects are continuous, they can 
only move the locations of the existing defects or, create 
or destroy defect-antidefect pairs. This is similar in spirit 
to the coarsening of defects observed in condensed mat- 
ter systems, though the actual dynamics of polarization 
defects is likely to be very different. 

In conclusion, we have drawn an analogy between the 
CMB polarization map and nematic liquid crystals, clas- 
sified the defects in the polarization map, and made a 
prediction for the scaling of the index of the CMB polar- 
ization map. Ongoing observations of the CMB will be 
able to find the exponent f3 and compare it to the theo- 
retical prediction of 1/2. Further statistical measures of 
the distribution of defects can provide additional novel 
probes of the CMB. 
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